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Abstract 

Values for the cross sections afj^ and crf^ t N are of crucial importance, once one 
wants to understand nuclear suppression of J/^(*S>') in a search for the quark gluon 
plasma. Quoted values are usually extracted from jN — ► J/*$>(^')N data via the 
vector dominance model (VDM) and are very small compared with expectations from 
QCD. The validity of QCD is questioned. When the VDM is extended to a multi- 
channel case, the 7iV data lead to a value 2.8 ± 0.12 mb < <T%£(y/a = 10 GeV) < 
4.1 ± 0.15 mb, which is a factor 3 larger than the VDM prediction. Errors are due 
to experimental input, and the corridor shows the theoretical uncertainty. We also 
derive from the data crf^ /af/^ N ~ 4, where we could not estimate the theoretical 
uncertainty. 
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1. Introduction 

The success of the vector dominance model (VDM) for the description of high-energy 
7P interactions, particularly for the production of light vector mesons ]T|| has led to a wide 
spread confidence that the VDM also accurately describes the production of heavy flavoured 
vector mesons. Using data on elastic photoproduction of charmonia, J/^ and the use 
of the VDM leads to a J/\I/-nucleon total cross section af Q ^ ~ 1.3 mb for a/s = 10 GeV and 
a fot N ' l a tJt ~ 0-8 ( see below), which are factors 2 to 4 smaller than predicted by QCD 0, ||. 

Reliable values for the charmonium interaction cross sections are of crucial importance, 
since J/\l/ and suppressions in heavy ion collisions are possible signatures for the quark- 
gluon plasma (QGP) formation. One needs to know the cross section af f in the absence 
of a QGP-environment in order to predict a nuclear suppression of J/^, which can be 
used as a base line in search for new physics. An analysis of available data for proton- 
nucleus collisions 0, [| in a simple model assuming instantaneous production of the J/^f |Q 
(which has a reasonable accuracy at medium energies ~ 10 GeV ||), leads to a value of 
a fJt ~ 6 mb |7|, which is substantially higher than the VDM prediction (see also ||). 

In this paper we argue that the discrepancy between the values for afj* from photopro- 
duction and hadroproduction data is to be sought in a too naive analysis of photoproduction 
data. 

Also a disagreement of the VDM prediction for the <pN cross section with the prediction 



of the additive quark model is claimed in || jTof. In this case more information about a 



<j>N 
tot 



is available than for charmonium. Facing such a problem with VDM in the strange sector, 



one reasonably concludes |10[ that predictions of the VDM for heavier flavours are even less 
trustable. 

The purpose of this paper is to clarify why the VDM fails for heavy flavours, and how 
it can be modified so that one can extract reliable values for heavy quarkonium interaction 
cross sections from photoproduction data. 

2. Predictions of the VDM 

The VDM establishes a relation between the elastic photoproduction 7^ — > VN cross 



section of a vector meson V and the total cross section a^J: of this meson on a nucleon [II], 



a 



tot J VDM 



l67ra em M v da(-yN -> VN) 



t=o 



1/2 



(1) 



3r e V(i + P^) 

Here py is the ratio of real to imaginary parts of the forward VN elastic scattering ampli- 
tude. In most cases, data are available only for the angle integrated cross section, which 
is related to the forward elastic photoproduction cross section by da^N — ► VN)/dt\ t =o = 
B^i N <j{"fN —>■ VN). The slope Bf { N of the elastic J/tyN differential cross section is de- 
termined mainly by the nucleon formfactor, since the radius of heavy quarkonium is small, 
(r 2 )^ <C (r 2 )N- Therefore, B^ N is only slightly larger than the half of the slope of elastic 



pp scattering. Our estimate B el 



6 GeV JTTJ] agrees well with experimental data \V1 



Refs. [14|, [T5J give a collection of experimental data, which we fit with the expression, 

2A 



a(jN -> VN) = a 



10 GeV 



(2) 



with parameters (T = (10.3±0.7) nb, X = 0.40±0.02. We can estimate = (7r/2)d[ln(a t *f )/dpn(s)] 
ttA/4 « 0.3. 

Applying the VDM relation (|l|) to the parameterization @ we extract a valueQ. 



( 



a ^)vDM = ( L2 4±0.13) mb x (^i/10 GeV) . 



(3) 



Using the VDM relation fll]) one can also extract the total cross section for the radial 
excitation ty'(2S). One finds for the ratio of the ^' to \& cross sections 



°tot 
a tot 



(4) 



with the ratio of the photoproduction cross sections R = a(jN — > ^/'N)/a( , ~fN — > tyN). 

gives 



For the experimental value R = 0.195 ± 0.03 |16[-]18[, eq. 



a 



tot 



°tot ) VDM 



0.8 ±0.2 . 



(5) 



This result contradicts dramatically the expectation based on QCD, 



"Actually, the exponent can be about 20% larger than A due to energy dependence of the slope Bf t N 



(6) 



This ratio is 7/3 in the harmonic oscillator model for J/ty and and equals 4 in a more 



realistic potential model [19 



3. Why does VDM fail for heavy flavours? 

The basic assumption of the VDM is that via quantum mechanical fluctuations the pho- 
ton converts into a vector meson V, which interacts with the proton elastically and thereby 
comes on the energy shell. Is the hadronic fluctuation of the 7 always a vector meson? 
In ref. |2(], 0, |2l| it is proposed that the photoproduction of J/\l/ should be considered in 
perturbative QCD as an interaction of a cc fluctuation (and not a J/$ fluctuation) of the 
photon with the target followed by projection of the produced cc wave packet on the J/ty 
wave function. Only if the wave function of the cc component of the photon is similar to 
that for the J/^f, this model is close to the conventional VDM. 

In order to see when and to what degree the wave functions of QQ fluctuations coincide 
with those of the vector mesons, considered as QQ bound states, we study the mean squared 
transverse radii. The wave function of the QQ component of the photon has the form M, 

% Q {r T ) oc K (m Q r T ) , (7) 

where K (x) is the modified Bessel function and r T is the transverse QQ separation. Its 
mean transverse size squared is calculated from (|7|) 

This value has to be compared with the mean transverse distance between the Q and Q in 
the vector meson, considered as the lowest state in an harmonic oscillator potential with 
frequency uj: 



(r 2 T )v = • (9) 

IllQ CU 
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Empirically one finds that ui = (My>(2S) ~ ^V(iS))/2 has a value of roughly 300 MeV 
independent of the flavour. 

The two sizes eqs. (§) and (|9|) depend differently on ttiq, they cannot coincide for all 
flavours, in particular, the discrepancy increases the heavier the quark mass is. 

We introduce the factor Nq, and for the transition from the VDM prediction to the 
real value 



°" = iV*^ J VDM (10) 

(and similarly for ty') and attempt to calculate these factors by two different methods in 
the next sections. 

4. Hadronic representation: a multi-channel approach 

The above observation that the cc fluctuation is quite different from the J/\I/ can be 
phrased formally in writing an expansion of the cc fluctuation in terms of a complete set of 
charmonium states 

|7> ce = ai|J/tf> + a 2 |tf , > + ... • (11) 

With this expansion of the cc fluctuation, photoproduction of and J/^f gains contributions 
from different mechanisms as is shown schematically in Fig. 2 for a two channel example: 
i) direct production as in the VDM 7 — > ty, — > and ii) indirect production, 7 — > 

^'iV — > ^iV, which contains the off-diagonal diffractive interaction # ^'iV. The two 
amplitudes which lead to the final state must be added before squaring. The non-diagonal 
process \I/ — > ^' is the correction to the VDM expression which we propose. 

The virtual photoproduction amplitude in the naive VDM reads, 



fvD M {l*N - *iV) = /(*JV - *iV) , (12) 



where the photon of virtuality Q 2 converts into with a probability amplitude C 1 ^ , the 
j/\[r propagates with a usual propagator and ffifN —>■ $?N) denotes the elastic scattering 
amplitude. The extension of eq. (|12|) to the multi-channel case is straightforward: 
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Figure V. Diagrams for the elastic photoproduction of the vector mesons, 
an example of two channels. Photoproduction of the ground state V (a), 
and of the radial excitation V (h). 



f( 7 *N - *N) = Nq,(Q 2 ) f VDM (j*N -> *N) , (13) 
where we have defined the correction factor Nq,(Q 2 ) 

N _ i | v f(%N -> *n) Ml + q 2 

iV*(Q)-l + ^ C i* fm _+ m Ml t + Qf (14) 
It follows from perturbative QCD (and is shown in the next section) that iV# oc 1/Q 2 for 
Q 2 — > oo. This property of color transparency implies for the expression ([14]) a sum rule, 

For heavy quarkonia the mass splitting is much smaller than the mass, therefore, we expect 
a strong cancellation in (|I1) even at low values of Q 2 . Indeed, the zero order term in an 



expansion of N<$(0) in the small parameter AM^/Mq, vanishes according to fllpp. The first 
nonzero term in the expansion is of the order of 

N * {0) ~ m^ 550 ' 2 ' (16) 



The observed strong cancellation between different channels in (|T3|) is a general property 
of heavy quarkonia, and we do not expect the factor N^(Q 2 ) to depend much on the set 
of orthogonal states chosen for the hadronic basis. In one case the function N^(Q 2 ) can 



be calculated exactly. If the operator for the elastic scattering amplitude of a cc pair on 
a nucleon is proportional to the transverse quark separation squared rf*, and harmonic 
oscillator wave functions are chosen for the hadronic basis, one has /(^JV — > ^iV) oc 
(\l/j|rf.|\l/o) — f° r 2 > 2. Then, according to the sum rule (|HJ), one finds the surprisingly 
simple result, 

w . (0 >) = , (17) 



which is independent of any matrix element. With N^(0) = 0.30, one arrives at a factor of 
3.3 correction to the naive VDM, 

, 0.4 



a*? = (4.3 ± 0.5) mb x 



10 GeV 



How good in the correction factor (fCTI)? We will give an independent calculation of 
Ny(Q 2 ) in the next sections. In fig. 2 we compare our predictions with the factor Nq,(Q 2 ) 
Eq. (|17]) and with N^(Q 2 ) = 1, which is the conventional VDM, and available data from 
the EMC experiment ||17|| . The comparison obviously supports the form (|17D , although the 
error bars are rather large. 

For the ^' production we have to introduce at least three intermediate states ^' and 
ty", which exhaust the sum in the hadronic oscillator basis and have 



x 



f(<y*p -> V'p) = f V DAl(j*P -> Vp) 

+ /(* f JV -> tf'N) M 2 , + Q 2 



(19) 



^ 2 C"r*' /(*'JV -> #W) M|. + Q 2 

where the squared bracket defines iVV (Q 2 ). In an oscillatory basis ffif'N — > ^'iV) / ffifN — ► 
\I/iV) = 7/3. Together with the color transparency sum rule (|15|) all the parameters are 
determined 



iMQ 2 ) = 

which expression leads to iVV(0) 



1 - 



2 M 2 , + Q 2 5 M 2 , + Q 2 



7 Ml, + Q 2 



(20) 



7 M 2 +Q 2 , 
0.064. Then the ratio of the corrected values for the 
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0.01 1 1 
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Q 2 (GeV 2 ) 

Figure 2: Data on the Q 2 -dependence of the exclusive muoproduction of 
J/^> J77|/ . The dashed curve is the expectation of the VDM, the solid 
curve includes the correction factor (|77[). 



ty'N to tyN total cross sections is 

= 3.75 . (21) 



hot iv^(UJ / o tot 



This ratio corresponds nearly exactly to the QCD prediction, namely, the ratio of radii 
squared, which calculated for a realistic potential |19[ gives a value 4. 



5. Quark representation 

In order to estimate the theoretical uncertainty of the correction factor (0) and N%> (0) 
evaluated using the harmonic oscillator basis, we calculate these factors by yet another 
approach. 



S 



The multi-channel equations (13) - (14) can be represented in the quark basis as, 

ftfN^*N) = (*\<T(r T )\'fJ (22) 

with the wave function for the photon |7* g ) in the quark basis. In the same representation, 
the expression for the VDM reads 

/( 7 *iV - *JV) = (*Kr T )|*)(*| 7c * a ) ■ (23) 
Then the correction factor takes the form 

NM > - (»K„)t)(»bsJ ' (24) 

where \l/ may stand for the wave functions of the J/\l/ or the The perturbative wave 
function of the cc fluctuation of the virtual photon reads in non-relativistic approximation, 
\j* 5 ) oc K (er T ) |22|, 0, where e 2 = m 2 + Q 2 /4. This perturbative QCD expression for |7* g ) 
neglects the interaction between c and c and therefore is questionable for low Q 2 . In this 
respect the hadronic representation is preferable, since effectively it takes into account all 
the non-perturbative effects. 

The mean squared cc separation in the photon is calculated to (r^) 7 * = 2/3(m 2 + Q 2 /4) 
and goes to zero for Q 2 — > oo. In this limit and using <r(rr) = ar 2 -,, the expression for Nq,(Q 2 ) 
becomes very simple, because ^(vt) can be replaced by ^(0) in the integrals involving 
K (r T ). Then, 

-> P^f (Q 2 - oo) . (25) 

This result is exact. We note, that (r^) 7 * goes to zero with the Q 2 — > oo and this leads to 
the sum rule fll5"D. If one uses this expression also at Q 2 = together with (r 2 -)^ = 2/um c 
eq. (H) one finds iV*(0) = 2u/m c , which agrees with (|17|) to lowest order in u/2m c . In the 
same limit 

AMQ 2 ) = (4k (06) 
iV*(Q 2 ) (r 2 )* • 1 ] 

The expression of N$,(Q 2 ) in (|25| ) is based on the approximation (r|,) 7 -C (r 2 )^, which 
is not well satisfied for charmonium. In order to test the accuracy, we give up this ap- 
proximation, and use another trial wave functions among them the harmonic oscillator 
ones. 



Be specific, we use a mixture of the harmonic oscillator (HO) and Coulomb (C) wave 
functions, 

= a(r\^ HO )+^(^c) (27) 
for IS* and 2S states. The Coulomb wave function is calculated with a running QCD 



coupling as is explained in |p3| |. We vary a {(3 = yl — a 2 ) and adjust the parameters of 



\^ho) and \^c) so that one has always the same (rf.)* as taken from . 

Then we obtain the following results: 
i) for the J/ty 

2.3 < A^ 1 < 3.3 , (28) 

where the upper limit corresponds to fll7D and the lower one is calculated from eq. (19bb) 
with a = in (|3). The case of pure Coulomb wave function (a = 0) lies in the interval. 



Thus, 

2.8 ± 0.12 < aff{Vs = 10 GeV) < (4.1 ± 0.15) mb (29) 
with the energy dependence 

*ZnV~s) = 4f (10 GeV) (j^) 0A (30) 



ii) The results for iVV(0) are very sensitive to the wave function (|27|) (because of the 
node), so that we prefer to give the result in the form ( f2~T|) without a corridor of theoretical 
uncertainty. 

Furthermore, we test the sensitivity of the results to the form of the transition operator 
cr(rr). We choose the form o"(r-r) = a[l — exp(— br^)]/b |[23[| , which gives the linear r\ 
dependence at small tt- We found a very small deviation of the order of few percent. 

6. Conclusion and discussions 

The conventional expressions of the VDM applied to the photoproduction data of char- 
monium production lead to inconsistencies both with experimental data (er*f values ex- 
tracted from hA data) as well as predictions from QCD (cr^ oc (r 2 )h)- A generalization of 
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the VDM to a multi-channel problem leads to a value af Q ^ about three times as large as the 
one predicted by VDM, and the discrepancies relating to af Q f are largely removed. 
Few more comments are in order. 

• Theoretical accuracy of the correction factor Nq,(Q 2 ) is estimated in the quark repre- 
sentation by varying the shape of the J/^J wave function and the form of the dipole 
cross section. We conclude that the uncertainty does not exceed 20%. 

• Another source of information about af Q f is production of J/\I / off nuclei, which give 
a slightly higher value than we predict. This is because the simplified analysis || 
does not take into account extension of the charmonium length of path in nuclear 
matter, due to the coherence length |J, which is about 2 fm for the kinematics of ||. 
Therefore a reduction of about 20% should be applied to Besides, dynamics of 
charmonium hadroproduction || is far more complicate than for photoproduction. 

In order to be free of the effects of the coherence and formation lengths one should 
produce a low energy charmonium off nuclei. A unique opportunity is annihilation 
pp — > \1> on a bound protons. 

• Using the same formulas for photoproduction of bottomia we arrive at the result, 
oUt = 8 \ a tot^) D ■ The cross section {o^j D as wen as the ratio r-j cannot be 
predicted, since no data on bottomium photoproduction is available yet. 

• We expect a weak rising energy- dependence of the correction factor N^(Q 2 ). Indeed, 
it is known from HERA data on the proton structure function F 2 (x, Q 2 ) that the 
smaller is the size of a photon fluctuation, the steeper is the energy dependence. This 
is why the J/\l/ photoproduction cross section growth so steeply compared to what 
is known from soft hadronic interactions. According to the QCD evolution equations 
this effect is due to fast grows of the gluon cloud of the cc fluctuation of the photon, 
which steeper than that of J . 

• One may wonder why VDM works so well for p mesons. Indeed, the same expression 
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leads to deviation of 40% from the VDM prediction. This estimate contains an 
additional uncertainty since the dipole approximation a{r T ) oc r\ is justified only at 
small while large distances are important for p mesons. A harmonic oscillator is 
also too rough an approximation in this case. 

• One should be cautious applying VDM to the photoproduction of J/\I/ on nucleons at 
low energy. If the lifetime of the J fluctuation in the photon is much shorter than 
the nucleon radius, t c ~ M|/2i? 7 <C R N , the photoproduction cross section should be 
much smaller than VDM predicts, since the fluctuation can interact only during its 
short lifetime. This effect is taken into account by the nucleon form factor (compare 
with |H). 
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